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Abstract. The concept of active Brownian particles is used to model a collective opinion formation process. 
It is assumed that individuals in community create a two-component communication field that influences 
the change of opinions of other persons and/or can induce their migration. The communication field is 
described by a reaction-diffusion equation, the opinion change of the individuals is given by a master 
equation, while the migration is described by a set of Langevin equations, coupled by the communication 
field. In the mean-field limit holding for fast communication we derive a critical population size, above 
which the community separates into a majority and a minority with opposite opinions. The existence of 
external support (e.g. from mass media) changes the ratio between minority and majority, until above 
a critical external support the supported subpopulation exists always as a majority. Spatial effects lead 
to two critical "social" temperatures, between which the community exists in a metastable state, thus 
fluctuations below a certain critical wave number may result in a spatial opinion separation. The range 
of metastability is particularly determined by a parameter characterizing the individual response to the 
communication field. In our discussion, we draw analogies to phase transitions in physical systems. 



PACS. 05.40.-a Fluctuation phenomena, random processes, noise, and Brownian motion 
organized systems - 87.23. Ge Dynamics of social systems 
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^ ■ 1 Introduction 
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In recent years, there has been a lot of interest in applica- 
tions of physical paradigms to a quantitative description 
of social and economic processes |^-[ll[| 



Methods of synergetics 
fl4| , IE], deterministic chaos 
els [20 " " 




stochastic processes 
and lattice gas mod- 
22| have been successfully applied for this purpose. 
The formation of public opinion p3|-p9| is among the 
challenging problems in social science, because it reveals a 
complex dynamics, which may depend on different inter- 
nal and external influences. We mention the influence of 
political leaders, the biasing effect of mass media, as well 
as individual features, such as persuasion or support for 
other opinions. 

A quantitative approach to the dynamics of opinion 
formation is given by the concept of social impact |2fj|,|23|], 
which is based on methods similar to the cellular automata 
approach |2^,^|. The social impact describes the force on 
an individual to keep or to change its current opinion. A 
short outline of this model is given in Sect. 2. The equi- 
librium statistical mechanics of the social impact model 
was formulated in pQj, while in 21,27,291 the occurrence 
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of phase transitions and bistability in the presence of a 
strong leader or an external impact have been analysed. 

Despite these extensive studies of the social impact 
model, there are several basic disadvantages of the con- 
cept. In particular, the social impact theory assumes, that 
the impact on an individual is updated with infinite veloc- 
ity, and no memory effects are considered. Further, there 
is no migration of the individuals, and any "spatial" dis- 
tribution of opinions refer to a "social" , but not to the 
physical space. 

In fact, the model of social impact has not been de- 
veloped to describe processes of opinion diffusion and mi- 
gration. In this paper, we present an alternative approach 
to the social impact model of collective opinion forma- 
tion, which tries to include these features. Our model is 
based on active Brownian particles, which interact via a 
communication field. This field considers the spatial dis- 
tribution of the individual opinions, further, it has a cer- 
tain life time, reflecting a collective memory effect and it 
can spread out in the community, modeling the transfer 
of information. 

Active Brownian particles ]3l]-|33| are Brownian parti- 
cles with the ability to take up energy from the environ- 
ment, to store it in an internal depot |m],|35| and to con- 
vert internal energy to perform different activities, such as 
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metabolism, motion, change of the environment, or signal- 
response behavior. As a specific action, the active Brown- 
ian particles (or active walkers, within a discrete approxi- 
mation) are able to generate a self-consistent field, which 
in turn influences their further movement and physical or 
chemical behavior. This non-linear feedback between the 
particles and the field generated by themselves results in 
an interactive structure formation process on the macro- 
scopic level. Hence, these models have been used to sim- 
ulate a broad variety of pattern formation processes in 
complex systems , ranging from physical to biological and 
social systems ]3rj-|4l|]. 

In Sect. 2, we specify the model of active Brownian 
particles for the formation of collective opinion structures. 
In Sect. 3, we discuss the limiting case of fast communica- 
tion between the individuals. Further, we investigate the 
influence of an external support and derive critical param- 
eters for the existence of subpopulations as majorities or 
minorities. In Sect. 4, we investigate spatial opinion struc- 
tures, and estimate critical wave numbers for the fluctua- 
tions, which lead to a spatial separation of the opinions. 
By deriving two different critical temperatures, we draw 
an analogy to the theory of phase transitions. 



2 Stochastic Model of Opinion Change and 
Migration 

Let us consider a 2-dimensional spatial system with the 
total area A, where a community of N individuals (mem- 
bers of a social group) exists. Each of them can share one 
of two opposite opinions on a given subject, denoted as 
9i = ±1; i = 1, N. Here, 9i is considered as an individ- 
ual parameter, representing an internal degree of freedom. 
Within a stochastic approach, the probability pi (9i , t) to 
find the individual i with the opinion 9i, changes in the 
course of time due to the following master equation: 

jUi (M) = $>W** )Pi $ , t) - 1* (M) X> $ 

6' 9'. 

(1) 

Here, means the transition rate to change the 

opinion 9i into one of the possible opinions 9[ during the 
next time step, with w(9i\9i) = 0. In the considered case, 
there are only two possibilities, either 9i = +1 — > 9\ = — 1, 
or 6i = — 1 — ► 9\ = +1. In the social impact theory |po| , p3|| , 
it is assumed that the change of opinions depends on the 
social impact, 7j, and a "social temperature", T [^l|,|27|]. 
A possible ansatz for the transition rate reads: 

w(9' i \9 i ) = r) exp{Ii/T}. (2) 

Here, r\ [1/s] defines the time scale of the transitions. T 
represents the erratic circumstances of the opinion change: 
in the limit T — > the opinion change is more determined 
by Ii, leading to deterministic transitions. As eq. (^) in- 
dicates, the likelyhood for changing the opinion is rather 
small, if Ii < 0. Hence, a negative social impact on indi- 
vidual i represents a condition for stability. To be specific, 



in the social impact theory, Ii may consist of three parts: 
I i = If + If+i r (3) 

If represents influences imposed on the individual by 
other members of the group, e.g. to change or to keep 
its opinion. If, on the other hand, is kind of a self-support 
for the own opinion, If < 0, and If x represents external 
influences, e.g. from government policy, mass media, etc. 
which may also support a certain opinion. 

Within a simplified approach of the social impact the- 
ory, every individual can be ascribed a single parameter, 
the "strength", Sj. Furthermore, a social distance dy is 
defined, which measures the distance between each two 
individuals in a social space [p0|,p3"|], which does not 
necessarily coincide with the physical space. It is assumed 
that the impact between two individuals decreases with 
the social distance in a non- linear manner. The above as- 
sumptions are included in the following ansatz pl|,p7[|: 

N 

Ii = -9i s A/ d ij - + e i i ( 4 ) 

e is the so-called self-support parameter, and n > is 
a model constant. The external influence, &i may be re- 
garded as a global preference towards one of the opinions. 
A negative social impact on individual i is obtained, (i) if 
most of the opinions in its social vicinity match its own 
opinion, or (ii) if the impact resulting from opposite opin- 
ions is at least not large enough to compensate its self- 
support, or (iii) if the external influences do not force the 
individual to change its opinion, regardless of self-support 
or the impact of the community. 

In the form outlined above, the concept of social im- 
pact has certain drawbacks: The social impact theory as- 
sumes that the impact on an individual is instantaneously 
updated, if some opinions are changed in the group (which 
basically means a communication with infinite velocity). 
Spatial effects in a physical space are not considered here, 
any "spatial" distribution of opinions refers to the social 
space. Moreover, the individuals are not allowed to move. 
Finally, no memory effects are considered in the social im- 
pact, the community is only affected by the current state 
of the opinion distribution, regardless of its history and 
past experience. 

In this paper, we want to modify the theory by includ- 
ing some important features of social systems: (i) the exis- 
tence of a memory, which reflects the past experience, (ii) 
an exchange of information in the community with a finite 
velocity, (iii) the influence of spatial distances between in- 
dividuals, (iv) the possibility of spatial migration for the 
individuals. It seems more realistic to us that individuals 
have the chance to migrate to places where their opinion is 
supported rather than change their opinion. And in most 
cases, individuals are not instantaneously affected by the 
opinions of others, especially if they are not in their close 
vicinity. 

As a basic element of our theory, a scalar spatio- 
temporal communication field hg(r,t) is used. Every in- 
dividual contributes permanently to this field with its 
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opinion 9i and with its personal strength Sj at its cur- 
rent spatial location r^. The information generated this 
way has a certain life time 1//3 [s] , further it can spread 
throughout the system by a diffusion-like process, where 
Dh [m 2 /s] represents the diffusion constant for informa- 
tion exchange. We have to take into account that there 
are two different opinions in the system, hence the com- 
munication field should also consist of two components, 
9 = {— 1, +1}, each representing one opinion. For simplic- 
ity, it is assumed that the information resulting from the 
different opinions has the same life time and the same way 
of spatial distribution; more complex cases can be consid- 
ered as well. 

The spatio-temporal change of the communication 
field can be summarized in the following equation: 



d_ 

df 



N 



hg(r, t) = Si 5 gfii S(r - r { 



i=l 



/3h e {r,t) + D h Ah e (r,t). 



(5) 



Here, 8g^ i is the Kronecker Delta indicating that the in- 
dividuals contribute only to the field component which 
matches their opinion 0j. S(r — r.j) means Dirac's Delta 
function used for continuous variables, which indicates 
that the individuals contribute to the field only at their 
current position, r,. We note that this equation is a 
stochastic partial differential equation with 



N 

£ 

i=l 



S(r — r.i) 



(0) 



being the microscopic density [B2| of the individuals 
changing their position due to Eq. (M). Hence, the changes 
of the communication field hg(r, t) are measured in units 
of a density of the personal strength Si. 

Instead of a social impact, the communication field 
hg(r,t) influences the individual i as follows: At a cer- 
tain location t^, the individual with opinion 6i = +1 
is affected by two kinds of information: the informa- 
tion resulting from individuals who share his/her opinion, 
£), and the information resulting from the oppo- 
nents hg=-i(ri,t). The diffusion constant Dh determines 
how fast he/she will receive any information, and the de- 
cay rate /3 determines, how long a generated information 
will exist. Dependent on the ZocaHnformation, the individ- 
ual has two opportunities to act: (i) it can change its opin- 
ion, (ii) it can migrate towards locations which provide a 
larger support of its current opinion. These opportunities 
are specified in the following. 

For the change of opinions, we can adopt the transition 
probability, Eq. (0) , by replacing the influence of the social 
impact Jj with the influence of the local communication 
field. A possible ansatz reads: 

w(9! i \9 i ) = V exp{[h e ,{r l: t) - h e (r l7 t)]/T} 
w(9i\9i) = (7) 

As in Eq. (||) , the probability to change opinion 9i is rather 
small, if the local field hg(ri,t), which is related to the 



support of opinion 9i , overcomes the local influence of the 
opposite opinion. This effect, however, is scaled again by 
the social temperature T, which is a measure for the ran- 
domness in social interaction. Note, that the social tem- 
perature is measured in units of the communication field. 

The movement of the individual located at space co- 
ordinate T{ may depend both on erratic circumstances 
and on the influence of the communication field. Within a 
stochastic approach, this movement can be described by 
the following overdamped Langevin equation: 



dt 



dh e (r,t) 



dr 



(8) 



In the last term of Eq. (j|) D n means the spatial diffusion 
coefficient of the individuals. The random influences on 
the movement are modeled by a stochastic force with a 
^-correlated time dependence, i.e. is the white noise 

with (&(t) &(*')} = <% <*(*-*!. 

The term h e (r,t) in Eq. (0) means an effective com- 
munication field which results from hg(r,t) as specified 
below. It follows that the overdamped Langevin Eq. (||) 
considers the response of the individual to the gradient 
of the field h e (r,t), where a, is the individual response 
parameter, weighting the importance of the information 
received. In the considered case, the effective communica- 
tion field h e (r, t) is a certain function of both components, 
h±i(r, t), of the communication field, see Eq. (||). One can 
consider different types of response, for example the fol- 
lowing: 

(i) The individuals try to move towards locations which 
provide the most support for their current opinion 6^. 
In this case, they only count on the information which 
matches their opinion, h e (r,t) — hg(r,t), and follow 
the local ascent of the field {pn > 0). 

(ii) The individuals try to move away from locations which 
provide any negative pressure on their current opinion 
9i . In this case, they count on the information resulting 
from opposite opinions (9'), h e (r,t) = h$<(r,t), and 
follow the local descent of the field (a^ < 0). 

(iii) The individuals try to move away from locations, if 
they are forced to change their current opinion 9i, 
but they can accept a vicinity of opposite opinions, 
as long as these are not dominating. In this case, they 
count on the information resulting from both support- 
ing and opposite opinions, and the local difference be- 
tween them is important: h e (r, t) — [h$(r, t) — hg'(r, t)] 
with on > 0. 

Additionally, the response parameter can also consider 
that the response occurs only, if the absolute value of 
the effective field is locally above a certain threshold h t hr- 
a>i = 0[\h e (r,i)\ — hthr], with 0[y] being the Heavyside 
function: = 1, if y > 0, otherwise = 0. We note that 
for the further discussions in Sect. 3 and 4, we assume 
h e (r,t) — hg(r,t) for the effective communication field 
(case i), while on = a is treated as a positive constant 
independent of i and h e (r,t). 

In order to summarize our model, we note the non- 
linear feedback between the individuals and the commu- 
nication field as shown in Fig. [I]. The individuals generate 
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the field, which in turn influences their further movement 
and their opinion change. In terms of synergetics, the field 
plays the role of an order parameter, which couples the in- 
dividual actions, and this way initiates spatial structures 
and coherent behavior within the social group. 




Fig. 1. Circular causation between the individuals with differ- 
ent opinions, C_i, C+i and the two-component communication 
field, he(r, t). 



The complete dynamics of the community can be for- 
mulated in terms of the canonical TV-particle distribution 
function 



P(i,L,t) =P(6 1 ,r 1 ,...,9 N ,r N ,t), 



(9) 



which gives the probability to find the N individuals with 
the opinions 9\, 9n in the vicinity of n, vn on the 
surface A at time t. Considering both opinion changes 
and movement of the individuals, the master equation for 
P{9 7 r, t) reads: 



d_ 

Of 



P(i,L,t) = [Hm')P(l',r,t)-w(9'\6)P(e,r,t) 
[V< (a Vihg{r, t) P(9,r, t)) - D n AP(r, 9, t) 



N 



i=l 



(10) 



The first line of the right-hand side of Eq. ( |l0| ) describes 
the "gain" and "loss" of individuals (with the coordinates 
ri,...,rjv) due to opinion changes, where w(9\9?) means 
any possible transition within the opinion distribution 9' 
which leads to the assumed distribution 9. The second 
line describes the change of the probability density due 
to the motion of the individuals on the surface. Eq. ( [To| ) 
together with eqs. (||), (0) forms a complete description of 
our system. 



3 The Case of Fast Communication 

3.1 Derivation of Mean Value Equations 

Let us first restrict to the case of very fast exchange of 
information in the system. Then, spatial inhomogenities 
are equalized immediately, hence, the communication field 
hg(r,t) can be approximated by a mean field hg{t): 



h e (t) = - I h.Ar. I x!r~ 



(11) 



where A means the system size. The equation for the mean 
field hg(t) results from eq. <M): 



dhgjt) 

dt 



-(3hg(t) + sng 



with Si = s and the mean density 



Ng_ 

A 



N 
A' 



(12) 



(13) 



where the number of individuals with a given opinion 9 
fulfils the condition 

^ N g = N +1 + N-i = N — const. (14) 



We note that in the mean-field approximation no spatial 
gradients in the communication field exist. Hence, there 
is no additional driving force for the individuals to move, 
as assumed in Eq. (|8|). Such a situation can be imagined 
for communities existing in very small systems with small 
distances between different groups. In particular, in such 
small communities also the assumption of a fast informa- 
tion exchange holds. Thus, in this section, we restrict our 
discussion to subpopulations with a certain opinion rather 
than to individuals at particular locations. 

Let p(Ng, t) denote the probability to find Ng individ- 
uals in the community which shares opinion 9. The master 
equation for p(N + i,t) explicitely reads: 

^P(N+I,t) = W(N +1 \N +1 - l)p(N +1 - l,t) 

+W(JV+i|JV+i + l)p(JV +1 + l,t) 

- P (N +1 ,t) |W(iV +1 + l|JV+i) 

+ W(N +1 -1\N +1 )]. (15) 

The transition rates W(M\N) appearing in Eq. ( p^| ) are 
assumed to be proportional to the probability to change 
a given opinion, eq. (Q) , and to the number of individuals 
which can change their opinion into the given direction: 

W(N +1 + l\N +1 )=N_ lV exp{(h +1 -h^ 1 )/T}, (16) 
W(N +1 - l\N +1 ) = N +1 r, exp {-(h +1 ~ h-i)/T}. 

The mean values for the number of individuals with a 
certain opinion can be derived from the master equation 

(PD 

(N e (t)) = Ngp(N e ,t), (17) 

where the summation is over all possible numbers of Ng 
which obey the condition eq. (|l4|). From eq. @, the de- 
terministic equation for the change of (Ng) can be derived 
in the first approximation as follows (4^] (see also [[[],§, (§): 



dt 



(Ng) = {W{Ng + l\Ng) - W(Ng - \\Ng)) (18) 



For iV+i, this equation reads explicitely: 



dt 



(N+i) = (N^ lV exp 



-N+i rj exp 



h+i(t) -h-i(t) 
T 

h+i{t) - 
T 



(19) 



Schweitzer/Holyst: Modelling Collective Opinion Formation 



5 



Introducing now the fraction of a subpopulation with opin- 
ion 9, xg — (Ng)/N, and using the standard approxima- 
tion to factorize eq. (|19|), we can write it as: 



x + i = (1 — x + i) r\ exp(a) — x+\ rj exp(— a), (20) 
a= \h +x (t) -h-tit)] IT. 

Via Ah(t) = h + \ —h—i, this equation is coupled with the 
equation 

Ah = -(3 Ah + sn(2x + i - l) (21) 



which results from eq. (12) for the two field components. 



3.2 Critical and Stable Subpopulation Sizes 

Within a quasistationary approximation, we can assume 
that the communication field relaxes faster than the dis- 
tribution of the opinions into a stationary state. Hence, 
with hg = 0, we find from Eq. (n2J) : 



L stat 
n +l 



13 x+1 



h stat 



K (x +1 - 



with k 



s n 
2s n 



1 



(22) 



Here, the parameter k includes the specific internal condi- 
tions within the community, such as the total population 
size, the social temperature, the individual strength of the 
opinions, or the life time of the information generated. In- 
serting a from eq. ( p2[ ) into eq. (22), a closed equation 
for xg is obtained, which can be integrated with respect 
to time (Fig. ||a). We find that, depending on k, different 
stationary values for the fraction of the subpopulations 
exist. For the critical value, k c = 2, the stationary state 
can be reached only asymptotically. Fig. ||(b) shows the 
stationary solutions, xg — 0, resulting from the equation 
for x+\\ 

(1 — x + \) exp [Kin] = x+\ exp [n (1 — x+i)] (23) 



For k < 2, x+i — 0.5 is the only stationary solution, 
which means a stable community where both opposite 
opinions have the same influence. However, for k > 2, 
the equal distribution of opinions becomes unstable, and 
a separation process towards a preferred opinion is ob- 
tained, where x±i = 0.5 plays the role of a separation 
line. We find now two stable solutions where both opin- 
ions coexist with different shares in the community, as 
shown in Fig. |[ Hence, each subpopulation can exist ei- 
ther as a majority or as a minority within the community. 
Which of these two possible situations is realized, depends 
in a deterministic approach on the initial fraction of the 
subpopulation. For initial values of x+i below the separa- 
trix, 0.5, the minority status will be most likely the stable 
situation, as Fig. ||(a) shows. 

The bifurcation occurs at k c — 2, where the former 
stable solution x +1 = 0.5 becomes unstable. From the 
condition k — 2 we can derive a critical population size, 



1.0 
0.8 
0.6 

-4— < 

X + 0.4 
0.2 
0.0 




0.0 



5.0 



10.0 



15.0 




+ 0.4 



iV c = (3AT/s, 



(24) 



Fig. 2. (a: top) Time dependence of the fraction x+i(t) (eq. 
pp| ) of the subpopulation with opinion +1 for different initial 
conditions and for three different values of k: 1.0 (solid line); 
2.0 (dot-dashed line), 3.0 (dashed line), (b: bottom) Station- 
ary solutions for x+\ (eq. ^) for different values of k. The 
bifurcation at the critical value k c = 2 is clearly visible. 



where for larger populations an equal fraction of opposite 
opinions is certainly unstable. If we consider e.g. a growing 
community with fast communication, then both contra- 
dicting opinions are balanced, as long as the population 
number is small. However, for N > N c , i.e. after a certain 
population growth, the community tends towards one of 
these opinions, thus necessarily separating into a majority 
and a minority. Which of these opinions would be dom- 
inating, depends on small fluctuations in the bifurcation 
point, and has to be investigated within a stochastic ap- 
proach. We note that eq. (pi]) for the critical population 
size can be also interpreted in terms of a critical social 
temperature, which leads to an opinion separation in the 
community. This will be discussed in more detail in Sect. 
4. 

From Fig. ||(b), we see further, that the stable coex- 
istence between majority and minority breaks down at a 
certain value of K, where almost the whole community 
shares the same opinion. From Eq. ( p3| ) it is easy to find 
that e.g. k « 4.7 yields x+i ~ {0.01; 0.99}, which means 
that about 99% of the community share either opinion +1 
or — 1. 
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3.3 Influence of External Support 

Now, we discuss the situation that the symmetry between 
the two opinions is broken due to external influences on 
the individuals. We may consider two similar cases: (i) 
the existence of a strong leader in the community, who 
possesses a strength s; which is much larger than the usual 
strength s of the other individuals, (ii) the existence of an 
external field, which may result from government policy, 
mass media, etc. which support a certain opinion with a 
strength s m . 

The additional influence s ex t '■= {sj/A,s m /A} mainly 
effects the communication field, eq. (|5|), due to an extra 
contribution, normalized by the system size A. 

If we assume an external support of opinion 9 = +1, 
the corresponding field equation in the mean-field limit 
(eq. |l2|) and the stationary solution (eq. ^2|) are changed 
as follows: 



h+i = -f3h + i(t) + snx + i + s ext 



h +1 = —X +1 + 



" = » ! x +i - 2 



~W' 



(25) 



Hence, in Eq. ( |23|) which determines the stationary solu- 
tions, the arguments are shifted by a certain value: 



(1 - x + i) exp 
x+i exp 



K X+i 



Sext 



\ S ex t 
K (1 — X+i 



(26) 



Fig. || shows how the critical and stable subpopulation 
sizes change for subcritical and supercritical values of k, 
dependent on the strength of the external support. 

For k < k c (Fig. |||a), we see that there is still only one 
stable solution, but with an increasing value of s ex t, the 
supported subpopulation exists as a majority. For k > n c 
(Fig. U b), we observe again two possible stable situations 
for the supported subpopulation, either a minority or a 
majority status. But, compared to Fig. ||(b), the symme- 
try between these possibilities is now broken due to the 
external support, which increases the region of initial con- 
ditions leading to a majority status. 

Interestingly, at a critical value of s ex t, the possibil- 
ity of a minority status completely vanishes. Hence, for a 
certain supercritical external support, the supported sub- 
population will grow towards a majority, regardless of its 
initial population size, with no chance for the opposite 
opinion to be established. This situation is quite often re- 
alized in communities with one strong political or religious 
leader ( "fundamentalistic dictatorships" ) , or in communi- 
ties driven by external forces, such as financial or military 
power ( "banana republics" ) . 

The value of the critical external support, s c ext , of 
course depends on k, which summarizes the internal situa- 
tion in terms of the social temperature, or the population 
size, etc. From Eq. (|26| ) we can derive the condition for 
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Fig. 3. Stable fraction of the subpopulation, as a func- 

tion of the strength s* = s ext /f3T of the external support, (a: 
top) k = 1, (b: bottom) k = 3. The dashed line in (b) repre- 
sents the separation line for the initial conditions, which lead 
either to a minority or a majority status of the subpopulation. 



which two of the three possible solutions coincide, thus 
determining the relation between s% xt and n as follows: 



^ = - In 

0T 2 



1 



1 - 



W 1 -! (27 » 



Fig. ^ shows how much external support is needed to 
paralyze a community with a given internal situation (n) 
by one ruling opinion. As one can see, the critical exter- 
nal support is an increasing function of the parameter k, 
meaning that it is more difficult to paralyse a society with 
strong interpersonal interactions. 

Let us conclude the discussion of the phase transition 
in the mean field limit, presented in this section. With re- 
spect to the social impact theory [p|pc|,pi],p3|,p9|, we note 
that phase transitions have been not considered there, 
since the focus was on other phenomena so far. On the 
other hand, for the case of two opinions, our results well 
correspond to those obtained by Weidlich and Haag in a 
model of collective opinion formation There, a mas- 
ter equation and appropriate utility potentials are used 
to analyse the mean-field dynamics of interacting popu- 
lations. Similar to our model, the approach used in 
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Fig. 4. Critical external support s* (eq. ^) as a function of 



leads to a phase transition and a corresponding bifurac- 
tion diagram. The main difference between both models 
is in the interpretation of the bifurcation parameter n. 
In our case, k results from other model parameters, e.g. 
from the mean "social strength" s which plays a role of a 
coupling constant between the opinions and the commu- 
nication field. In the model of Weidlich and Haag, on the 
other hand, this parameter is interpreted as a derivative of 
utility potentials. Our analytic result for the critical exter- 
nal support, Eq. (p7|), is in qualitative agreement with the 
stability analysis and the computer simulations presented 
in §§. 



4 Critical Conditions for Spatial Opinion 
Separation 

In the previous section, the existence of critical parame- 
ters, such as k c or s% xt , has been proved for a community 
with fast communication, where no inhomogenitics in the 
communication field can exist. In the more realistic case, 
however, we have finite diffusion coefficients for the in- 
formation, and the mean-field approximation, eq. (|l2|), is 
no longer valid. Instead of focussing on the subpopulation 
sizes, we now need to consider the spatial distribution of 
individuals with opposite opinions. 

Starting with the canonical N-particle distribution 
function, P(0,r,t), eq. JlCj), the spatio-temporal density 
of individuals with opinion 6 can be obtained as follows: 



N 



ng(r,t) = J^2$e,8 z 5(r-r 



xP(6 1 ,r 1 ...,9 N ,r N ,t) dr l ...dr N (28) 

Integrating eq. according to eq. ( |28| ) and neglecting 
higher order correlations, we obtain the following reaction- 
diffusion equation for ng(r, t) 



d_ 

di 



ng(r,t) 



ng(r,t)aVhg(r,t) + D n Ang(r,t) 



- J2 [w(6'\6)n e (r,t) +w(6\6')ne'(r,t) 

(29) 

with the transition rates obtained from eq. (0) : 

w(6'\6) = r] exp{[he>(r, t) - h e (r,t)]/T} (30) 
w(6\6) = 

With 9 = {+1,-1}, Eq. ( p9| ) is a set of two reaction- 
diffusion equations, coupled both via ng(r, t) and hg(r, t). 
Inserting the densities ng(r, t) and neglecting any external 
support, eq. (|b|) for the spatial communication field can be 
transformed into the linear deterministic equation: 

^-h e {r,t) = sn e (r,t) - fth e (r,t) + D h Ah s (r,t) (31) 
at 

The solutions for the spatio-temporal distributions of in- 
dividuals and opinions are now determined by the four 
coupled equations, Eq. ( |29| ) and Eq. (|3l|). For our fur- 
ther discussion, we assume again that the spatio-temporal 
communication field relaxes faster than the related distri- 
bution of individuals into a quasi-stationary equilibrium. 
The field hg(r, t) should still depend on time and space co- 
ordinates, but, due to the fast relaxation, there is a fixed 
relation to the spatio-temporal distribution of individuals. 
Further, we neglect the independent diffusion of informa- 
tion, assuming that the spreading of opinions is due to the 
migration of the individuals. From Eq. (|3l]), we find with 
h e (r,t) = and D h = 0: 



he(r,t) = -n e (r,t) 



(32) 



which can now be inserted into Eq. (j29|), thus reducing 
the set of coupled equations to two equations. 

The homogeneous solution for ng(r,t) is given by the 
mean densities: 



n g = (ne{r,t)) = - 



(33) 



Under certain conditions however, the homogeneous state 
becomes unstable and a spatial separation of opinions oc- 
curs. In order to investigate these critical conditions, we 
allow small fluctuations around the homogeneous state fig: 



ng(r,t) — n e + 5n g 



Sng 



no 



< 1 



(34) 



Inserting Eq. (B4|) into Eq. (E9T), a linearization gives: 



dSng 



D„ 



2ft 



ftT 

(35) 
(36) 

we find from Eq. ( |35| ) the dispersion relation A(fc) for small 
inhomogeneous fluctuations with wave vector k. This re- 
lation yields two solutions: 

Ai(fe) = -k 2 C + 2B ; \ 2 {k)=_~k 2 C 
n; C = D n - asfl 



With the ansatz 



Sn ~ exp (At + ikr) 



R - VSn 

ftT 



2ft 



(37) 
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For homogeneous fluctuations we obtain from Eq. (p 



The second critical temperature T| reads as follows: 



Ai 



2 r\ s n 

(3T 



2 77; A 2 = for k = 



(38) 



= _a_ sn = a_ 
2 2\i [3 2jji 1 



(44) 



which means that the homogeneous system is marginally 
stable as long as Ai < 0, or s h/f3T < 1. This result agrees 
with the condition k < 2 obtained from the previous mean 
field investigations in Sect. 3. The condition k = 2 or 
B = 0, respectively, defines a critical social temperature 



rpc _ STL 



(39) 



For temperatures T < Tf, the homogeneous state 
ng(r,t) — h/2, where individuals of both opinions are 
equally distributed, becomes unstable and the spatial sep- 
aration process occurs. This is in direct analogy to the 
phase transition obtained from the Ising model of a fer- 
romagnet. Here, the state with n < 2 or T > Tf, re- 
spectively, corresponds to the paramagnetic or disordered 
phase, while the state with n > 2 or T < Tf , respectively, 
corresponds to ferromagnetic ordered phase. 

The conditions of Eq. ( |38| ) denote a homogeneous sta- 
bility condition. To obtain stability against inhomoge- 
neous fluctuations of wave vector k, the two conditions 
Ai(fe) < and A 2 (fc) < have to be satisfied. 

Taking into account the critical temperature Tf , Eq. 
(|39|), we can rewrite these conditions, Eq. (|37|), as follows: 



D n - D[ 



k 2 [D„ 



— 1 J > 

D C S\ > 



(40) 



D'„ 



a sn 

2 T 



Hence, the condition 



D n > m 



(41) 



(42) 



fe 2 /i(r-T^ 



k 2 ^(r-T^ 



> 

> 



The occurence of two critical social temperatures Tf , T 2 
allows a more detailed discussion of the stability condi- 
tions. Therefore, we have to consider two separate cases of 
Eq. @: (1) Tf > Tf and (2) Tf < Tf , which correspond 
either to the condition a < 2/i, or a > 2/i, respectively. 

In the first case, Tf > T%, we can discuss three ranges 
of the temperature T: 

(i) For T > T{ both eigenvalues Ai(fe) and A 2 (fc), Eq. 
(pTj), arc nonpositive for all wave vectors k and the 
homogenous solution n/2 is completely stable. 

iii) For Tf > T > T 2 C the eigenvalue A 2 (fc) is still non- 
positive for all values of fc, but the eigenvalue Ai(fe) 
is negative only for wave vectors that are larger than 
some critical value k 2 > k 2 :. 



2 2 77 Tf — T 



k = 



fiT T — T 2 C 



(45) 



Here, a critical diffusion coefficient for the individuals 
appears, which results from the condition C = 0: 



denotes a second stability condition. In order to explain 
its meaning, let us consider that the diffusion coefficient 
of the individuals, D ni may be a function of the social 
temperature, T. This sounds reasonable since the social 
temperature is a measure of randomness in social inter- 
action, and an increase of such a randomness lead to an 
increase of a random spatial migration. The simplest re- 
lation for a function D n (T) is the linear one, D n = /iT. 
By assuming this, we may rewrite Eq. (fie] ) using a sec- 
ond critical temperature, T% instead of a critical diffusion 
coefficient WL: 



(43) 



This means that, in the given range of temperatures, 
the homogeneous solution n/2 is metastable in an in- 
finite system, because it is stable only against fluc- 
tuations with large wave numbers, i.e. against small- 
scale fluctuations. Large-scale fluctuations destroy the 
homogeneous state and result in a spatial separation 
process, i.e. instead of a homogenous distribution of 
opinions, individuals with the same opinion form sep- 
arated spatial domains which coexist. The range of the 
metastable region is especially determined by the value 
of a < 2/i, which defines the difference between Tf and 
T 2 C . 

(iii) For T < T| both eigenvalues Ai(fe) and A 2 (fc) are pos- 
itive for all wave vectors k (except k = 0, for which 
A 2 = yields), which means that the homogeneous so- 
lution n/2 is completely unstable. On the other hand all 
systems with spatial dimension L < 2ix /k c are stable 
in this temperature region. 

For case (2), Tf < T 2 , which corresponds to a > 2u, al- 
ready small inhomogeneous fluctuations result in an in- 
stability of the homogenous state for T < T 2 , i.e we 
have a direct transition from the completely stable to 
the completely unstable regime at the critical tempera- 
ture T = T 2 C . 

That means the second critical temperature T 2 marks 
the transition into complete instability. The metastable 
region, which exists for a < 2u, is bound by the two 
critical social temperatures, Tf and Tf. This allows us 
again to draw an analogy to the theory of phase transi- 
tions . It is well known from phase diagrams that the 
density-dependent coexistence curve Tf (n) divides stable 
and metastable regions, therefore we can name the critical 
temperature Tf , Eq. (|39|), as the coexistence temperature, 
which marks the transition into the metastable regime. On 
the other hand, the metastable region is separated from 
the completely unstable region by a second curve T 2 (n), 



Schweitzer/Holyst: Modelling Collective Opinion Formation 



9 



known as the spinodal curve, which defines the region of 
spinodal decomposition. Hence, we can identify the second 
critical temperature Tj, Eq. (|44|), as the instability tem- 
perature. 

We note that similar investigations of the critical sys- 
tem behavior can be performed by discussing the depen- 
dence of the stability conditions on the "social strength" 
s or on the total population number N — Ah. These in- 
vestigations allow the calculation of a phase diagram for 
the opinion change in the model discussed, where we can 
derive critical population densities for the spatial opinion 
separation within the community. 

5 Conclusions 

We have discussed a simple model of collective opinion 
formation, based on active Brownian particles, which rep- 
resent the individuals. Every individual shares one of two 
opposite opinions and indirectly interacts with its neigh- 
bours due to a communication field, which contains the 
information about the spatial distribution of the different 
opinions. This two-component field has a certain lifetime, 
which models memory effects. Furthermore, it can spread 
out in the community, which describes the diffusion of in- 
formation. This way, every individual locally receives in- 
formation about the opinion distribution, which affects its 
further actions: (i) the individual can keep or change its 
current opinion, or (ii) it can stay or migrate towards re- 
gions where its current opinion is supported. Both actions 
depend (a) on a social temperature, which describes the 
stochastic influences, and (b) on the local strength of the 
communication field, which expresses the deterministic in- 
fluences of the decision of an individual. 

For supercritical conditions within the community 
(e.g. supercritical population size, or supercritical external 
pressure, or low temperature etc.), the non-linear feedback 
between the individuals and the communication field, cre- 
ated by themselves, results in a process of spatial opinion 
separation. In this case, the individuals either change their 
opinion to match the conditions in their neighbourhood, 
or they keep their opinion, but migrate into regions which 
support this opinion. 

In this paper, we have studied the critical conditions, 
which may lead to this separation process. In the spatially 
homogeneous case, which holds either for small communi- 
ties or for an information exchange with infinite velocity, 
the communication field can be described in a mean-field 
approximation. 

For this case, we derived a critical population size, N° 
(which is related to a critical social temperature, T-f). For 
N < N c , there is a stable balance where both opinions are 
shared by an equal number of individuals. For N > N c , 
however, one of these opinions becomes preferred, hence, 
majorities and minorities appear in the community. Fur- 
ther, we have shown how these majorities change if we 
consider an external support for one of the opinions. We 
found, that beyond some critical support, the supported 
subpopulation must always exist as a majority, since the 
possibility of its minority status simply disappears. 



As a second case, we have investigated a spatially inho- 
mogeneous communication field, which is locally coupled 
to the distribution of the individuals. This coupling is due 
to an adiabatically fast relaxation of the communication 
field into a quasistationary equilibrium. 

Using this adiabatic approximation, we were able to 
derive critical conditions for a spatial separation of opin- 
ions. We found that above the critical population size 
(or for T < Tf), the community could be described as 
a metastable system, which expresses stability against 
small-scale pertubations. The region of metastability is 
bound by a second critical temperature, T 2 C , which de- 
scribes the transition into instability, where every pertu- 
bation results in an immediate separation. Further, we 
obtained that the range of metastability is particularly 
determined by the parameter a, which characterizes how 
strong an individual responses to the information received 
from the communication field. 

Finally, we would like to note that our model of col- 
lective opinion formation only sketches some basic fea- 
tures of structure formation in social systems. There is no 
doubt, that in real human societies a more complex be- 
havior among the individuals occurs, and that decision 
making and opinion formation may depend on numer- 
ous influences beyond a quantitative description. In this 
paper, we restricted ourselves to a simplified dynamical 
approach, which purposely stretches some analogies be- 
tween physical and social systems. The results, however, 
display similarities to phenomena observed in social sys- 
tems and allow an interpretation within such a context. 
So, our model may give rise to further investigations in 
the field of quantitative sociology. 
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